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This  report  describes  the  experimental  investigation  of  the  possibility 
of  real-time  Image  compensation  of  the  solar  photosphere  using  the  21- 
channel  Real-Time  Atmospheric  Compensation  (RTAC)  system  developed  at 
Itek. 

Observations  were  made  using  the  76-cm  Vacuum  Tower  Telescope  at 
Sacramento  Peak  Observatory  in  June  1979  and  March  1980.  During  the 
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second  period,  observations  were  also  made  of  single  and  double  stars  to 
obtain  data  on  the  isoplanatic  patch.  The  compensated  image  was  recorded 
on  videotape  using  an  IS1T  camera.  In  addition,  48  channels  of  real-time 
RTAC  data  (including  deformable  mirror  drive  signals)  were  digitized  and 
recorded  on  9-track  tape. 

With  stars,  the  real-time  compensation  system  operated  as  expected, 
sharpening  the  image  and  producing  a  significant  Increase  in  peak  Intensity 
Observations  of  double  stars  indicated  that  the  isoplanatic  patch  radius 
at  night  was  between  2  and  4  arcseconds  near  the  zenith. 

When  observing  the  sun,  it  was  found  that  the  image  stabilization  system 
worked  efficiently  even  with  bad  turbulence.  However,  the  expected  improve 
ment  in  image  sharpness  due  to  the  deformable  mirror  was  obtained  only  on 
some  occasions  and  could  not  generally  be  predicted,  in  spite  of  the  fact 
that  the  RTAC  appeared  to  operate  normally  with  all  feedback  loops  closed 
and  small  residual  errors. 

These  results  appear  to  be  due  to  atmospheric  parameters,  primarily  the 
small  size  of  the  isoplanatic  patch  and  the  small  values  of  turbulence 
coherence  length  during  the  day. 

These  two  factors  reduce  both  the  compensated  field  size  and  the 
compensated  aperture  size  (for  a  fixed  number  of  actuators)  resulting  in 
a  small  number  of  compensated  Image  elements. 

This  experiment  has  confirmed  the  feasibility  of  real-time  compensation 
of  solar  images.  However,  much  work  remains  to  be  done  to  refine  the 
technique.  Futher  operational  systems  should  employ  larger  apertures 
(necessitating  a  much  larger  number  of  actuators)  in  order  to  maximize 
the  number  of  resolved  image  elements  within  the  compensated  field. 
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1. 


OBJECTIVES 


The  prime  objective  of  this  program  is  to  make  an  experimental 
investigation  of  the  feasibility  of  real-time  compensation  of  the 
solar  photosphere,  with  the  goal  of  obtaining  high- resolution  images 
showing  the  structure  and  evolution  of  the  solar  granulation  pattern. 
Such  data  are  of  great  potential  value  in  advancing  our  understanding 
of  the  physical  processes  taking  place  in  the  sun.  An  important 
practical  consequence  will  be  an  improvement  in  our  capability  of 
predicting  magnetic  storms  and  radiation  outbursts  that  have  a  con¬ 
siderable  effect  on  military  communications  and  sensing  systems. 

The  best  currently  available  images  of  solar  granulation  obtained 
with  ground-based  telescopes  are  limited  by  atmospheric  turbulence 
to  a  resolution  of  1  to  2  arc-sec  with  occasional  frames  showing 
sub-arc-sec  detail.  Speckle  interferometric  measurements  made  with 
the  76-cm  Vacuum  Tower  Telescope  at  Sacramento  Peak  Observatory 
indicate  the  presence  of  detail  in  the  solar  granulation  up  to  at 
least  2.5  cycles/arc- sec .  An  initial  goal  was  therefore  to  demon¬ 
strate  consistent  imaging  of  0.4  arc-sec  solar  detail,  a  goal  that 
should  be  achievable  using  a  modified  version  of  Itek's  21-channel 
Real-Time  Atmospheric  Compensation  (RTAC)  system. 

The  RTAC  employs  a  deformable  mirror  with  21  independent  actuators, 
driven  in  the  closed-loop  mode  by  a  wavefront  sensor  which  operates 
with  white-light  reference  sources.  The  response  time  of  the  RTAC 
feedback  loop  is  about  1  millisecond,  which  is  fast  enough  to  com¬ 
pensate  atmospheric  turbulence  in  ground-based  telescopes.  For  the 
normal  range  of  atmospheric  turbulence  strength  at  good  sites 
(r0  *  6  to  12  cm),  the  RTAC  is  capable  of  compensating  telescope 
apertures  of  30  to  60  cm.  A  full  description  of  the  RTAC  is  given 
in  the  Reference  1. 

The  RTAC  was  originally  designed  to  operate  with  near-point  reference 
sources,  such  as  stars  and  other  smal 1 -diameter  astronomical  objects. 
Recent  laboratory  experiments  have  shown  the  feasibility  of  modifying 
the  wavefront  sensor  so  that  very  bright  objects  (such  as  the  sun) 
containing  detail  in  the  sub- arc- second  range  can  be  used  as  the 
source.  The  present  experiment  is  therefore  planned  to  validate  this 
application  of  the  RTAC  using  small  sections  of  the  solar  disk  as  the 
wavefront  reference. 

If  successful,  the  experiment  should  yield  improved  solar  images  over 
extended  times  of  several  hours.  It  should  also  allow  the  use  of 
high  spatial  resolution  spectroscopy  by  isolating  sections  of  the 
solar  granulation  that  are  normally  blurred  by  atmospheric  seeing. 

The  impact  of  this  technique  on  the  future  of  ground-based  solar 
observations  is  considerable,  as  it  will  give  a  new  lease  of  life  to 
the  many  existing  solar  observatories,  and  allow  the  collection  of 
high  resolution  data  otherwise  obtainable  only  from  space. 
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The  secondary  objective  of  the  program  is  to  obtain  measurements  of 
the  size  of  the  isoplanatic  patch  by  observing  binary  stars  of 
varying  separations.  To  perform  this  test,  the  RTAC  wavefront  sensor 
is  modified  by  the  addition  of  a  field  stop  to  isolate  the  light  from 
one  star  of  the  pair,  so  that  the  wavefront  measurement  and  resulting 
compensation  are  valid  for  only  one  line  of  sight.  By  recording  the 
images  of  both  stars,  and  determining  the  image  quality  of  the  second 
star  as  a  function  of  the  angular  separation,  the  size  of  the  iso¬ 
planatic  patch  can  be  measured. 

In  addition  to  the  compensated  images  recorded  on  videotape,  a  con¬ 
siderable  data  base  of  real-time  outputs  from  the  RTAC,  such  as 
deformable  mirror  drive  signals,  will  also  be  obtained. 
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EQUIPMENT 


2 . 

The  Real-Time  Atmospheric  Compensation  (RTAC)  system  modified  for 
solar  operation  was  installed  on  a  12-foot  optical  bench  on  the 
operating  platform  of  the  76-cm  Vacuum  Tower  Telescope  at  Sacramento 
Peak  Observatory  in  June  1979.  A  schematic  of  the  system  is  shown 
in  Figure  1.  The  prime  solar  image,  which  is  approximately  0.5 
meter  in  diameter,  was  focused  via  the  exit  port  and  fold  mirror 
at  the  end  of  the  optical  bench.  Most  of  the  sunlight  was  excluded 
by  a  field  stop  at  this  point,  only  a  small  area  of  about  20  x  20 
arcseconds  being  used  at  any  one  time.  The  area  viewed  was 
selected  by  driving  the  turret  mirrors  in  the  normal  manner.  The 
light  from  the  selected  image  area  was  collimated  to  form  a  beam 
exactly  1.00  inch  (25.4  mm)  in  diameter,  which  was  then  reflected 
from  the  tip-tilt  (guider)  and  deformable  mirrors  in  the  RTAC.  An 
optical  zoom  system  on  the  bench  varied  the  effective  focal  length 
of  the  collimator  to  allow  compensation  of  primary  apertures  between 
30  and  60  cm.  The  deformable  mirror  had  five  actuators  across  its 
diameter  so  that  the  effective  actuator  spacing  could  be  varied  in 
this  way  from  6  cm  to  12  cm  to  allow  for  difference  atmospheric 
turbulence  coherence  lengths. 

The  RTAC  system  was  modified  for  solar  operation  by  placing  a  field 
stop  in  the  wavefront  sensor  to  define  small  areas  of  the  solar 
surface,  adjustable  from  about  2  to  8  arcseconds  across,  for  use  as 
the  wavefront  reference  source.  The  wavefront  sensor  could  be 
adjusted  to  lock  on  to  either  a  bright  object  in  a  dark  field  or  a 
dark  object  in  a  bright  field. 

The  deformable  mirror  employed  in  the  RTAC  was  a  monolithic  piezo¬ 
electric  mirror  (MPM)  with  21  independently  controlled  zones.  The 
bandwidth  of  the  servo  control  loops  could  be  set  at  50,  200  or 
400  Hz.  The  two-axis  tip-tilt  mirror  used  for  overall  tilt  (or 
image  motion)  removal  had  a  bandwidth  of  150  Hz. 

To  provide  a  real-time  indication  of  RTAC  system  performance,  the 
residual  error  signals  in  all  21  control  loops  were  squared  and 
summed  to  form  a  composite  overall  error  voltage  that  was  displayed 
on  a  meter. 

The  imaging  system  consisted  of  a  525-line  TV  camera  with  relay 
optics  to  give  a  f ield-of- view  of  18  x  14  arcseconds  at  0.3  meter 
telescope  aperture,  or  9  x  7.6  arcseconds  at  0.6  meter  aperture. 
Composite  video  signals  were  recorded  on  analog  videotape  using  a 
video  cassette  recorder.  The  v'ideo  image  format  and  scale  factors 
for  the  range  of  field  stops  used  is  shown  in  Figure  2. 

A  real-time  digital  data  recording  system  was  attached  to  the  RTAC 
to  allow  up  to  48  channels  of  data  at  1  KHz  sampling  rate  to  be 
recorded  on  9-track  digital  tape.  Approximately  four  minutes  of 
continuous  observation  could  be  stored  on  each  2400  foot  tape. 

The  signals  recorded  are  listed  in  Table  1  and  included  the  drive 
signals  for  the  tip-tilt  mirror  and  MPM  actuators,  residual  wave- 
front  error,  and  wavefront  sensor  light  intensities. 
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MPM  -  Monolithic  piezoelectric  mirror 
TT  -  Tip- tilt  corrector 
WFS  -  Wavefront  sensor 
VC  -  Video  camera 
VM  -  Video  monitor 


Fig.  1  Solar  telescope  image  compensation  system 
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TABLE  1 


RTAC  DATA  CHANNEL  ALLOCATION 


Channel 


Number 

Description 

Signal 

1 

Tip/tilt  loop  status 

Ov-on,  +5v-off 

MPM  loop  status 

Ov-on,  +5v-off 

3 

X-tilt  drive 

*  lOv 

4 

Y-tilt  drive 

*  lOv 

5 

Command- start  to  record  data 

6 

MPM  drive  node  no.  5 

*  lOv 

7 

MPM  drive  node  no.  6 

±10v 

8 

MPM  drive  node  no.  7 

*10v 

9 

MPM  drive  node  no.  8 

*  lOv 

10 

MPM  drive  node  no.  9 

*  lOv 

11 

MPM  drive  node  no.  10 

*  lOv 

12 

MPM  drive  node  no.  11 

*  lOv 

13 

MPM  drive  node  no.  12 

*10v 

14 

MPM  drive  node  no.  13 

*10v 

15 

Seeing  monitor  output 

0-10  volts 

16 

Residual  error  output 

0-7  volts 

17 

MPM  drive  node  no.  14 

*10v 

18 

MPM  drive  node  no.  15 

*  lOv 

19 

MPM  drive  node  no.  16 

1 1 0  v 

20 

MPM  drive  node  no.  17 

1  lOv 

21 

MPM  drive  node  no.  18 

*  lOv 

22 

MPM  drive  node  no.  19 

*10v 

23 

MPM  drive  node  no.  20 

4  lOv 
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TABLE  1  (Continued) 


Chanel 


Number 

Description 

Signal 

24 

MPM  drive  node  no 

.  21 

lOv 

25 

MPM  drive  node  no 

.  1 

lOv 

26 

MPM  drive  node  no 

.  2 

lOv 

27 

MPM  drive  node  no 

.  3 

lOv 

28 

MPM  drive  node  no 

.  4 

lOv 

29 

Command  -  use  RTAC  tracking 

signals 

30 

X-tilt 

phase  error 

31 

Y-tilt 

phase  error 

32 

Light 

intensity, 

PMT 

XI 

33 

Light 

intensity, 

PMT 

X2 

34 

Light 

intensity, 

PMT 

X3 

35 

Light 

intensity. 

PMT 

X4 

36 

Light 

intensity , 

PMT 

X5 

37 

Light 

intensity. 

PMT 

X6 

38 

Light 

intensity , 

PMT 

X7 

39 

Light 

intensity , 

PMT 

X8 

40 

Light 

intensity. 

PMT 

X9 

41 

Light 

intensity , 

PMT 

X10 

42 

Light 

intensity , 

PMT 

Xll 

43 

Light 

intensity, 

PMT 

X12 

44 

Light 

intensity , 

PMT 

X13 

45 

Light 

intensity , 

PMT 

X14 

46 

Light 

intens i ty , 

PMT 

X15 

47 

Light 

intensity. 

PMT 

X16 

48 


Phase  error,  X4  Phase  detector  output 


The  relation  of  the  deformable  mirror  actuator  zones  to  the  detector 
locations  in  the  RTAC  is  shown  in  Figure  3. 

The  scale  factors  for  the  analog  and  digital  outputs  are  listed  in 
Table  2. 
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PHOTOMULTIPLIER 
DETECTOR  LOCATIONS 


DEFORMABLE 
MIRROR  ACTUATOR 
LOCATIONS 


FIGURE  3.  ACTUATOR  AND  DETECTOR  LOCATIONS 


TABLE  2 

RTAC  DATA  TAPE  CALIBRATION 


3. 


OBSERVATIONS 


Observations  were  made  using  the  RTAC  with  the  76-cm  Vacuum  Tower 
Telescope  at  Sacramento  Peak  National  Observatory  during  two  periods, 
18  through  28  June  1979  and  21  March  through  7  April  1980.  The 
telescope  time  that  was  allocated  enabled  observations  of  the  sun 
during  the  day  and  bright  stars  at  night. 

3.1  First  Observing  Period 

3.1.1  Stellar  Images 

It  was  found  that  the  RTAC  produced  a  significant  im¬ 
provement  in  stellar  images,  stabilizing  the  motion  of  the  seeing 
disk  and  reducing  its  size  in  many  cases  to  close  to  the  diffraction 
limit  of  the  telescope  (0.2  arc  seconds).  Compensated  and  uncom¬ 
pensated  images  of  the  Star  Arcturus  are  shown  in  Figure  4. 

3.1.2  Solar  Images 

When  operating  on  the  sun,  the  objects  used  as  reference 
sources  were  pores  and  small  sunspots  a  few  arcseconds  in  diameter. 
Attempts  to  lock  the  RTAC  on  to  the  solar  granulation  were  not 
success  ful . 

The  RTAC  achieved  image-motion  stabilization  using  the  tip-tilt 
mirror  without  any  difficulty  on  almost  all  occasions,  even  in  bad 
turbulence.  When  compensating  image  motion,  the  reference  object 
in  the  center  of  the  field  of  view  was  well  stabilized,  but  residual 
image  motion  due  to  isoplanatism  was  usually  observed,  increasing  as 
a  function  of  radius  from  the  center. 

The  result  of  operating  the  deformable  mirror  on  solar  images  was 
somewhat  ambiguous.  Again,  it  was  possible  to  close  the  servo  loops 
(21  in  all,  one  to  each  actuator  of  the  deformable  mirror)  on  most 
occasions,  even  in  bad  turbulence.  The  residual  wavefront  error,  as 
indicated  by  the  phase  error  variance  at  the  output  of  the  wavefront 
sensor,  was  usually  reduced  to  almost  zero,  showing  that  the  RTAC 
was  in  fact  efficiently  compensating  the  wavefront  errors  being 
measured.  However,  a  noticeable  improvement  in  the  resolution  of 
the  optical  image  was  apparent  only  on  rare  occasions.  The  most 
likely  reason  for  the  lack  of  improvement  in  image  detail  was  that 
the  isoplanatic  patch  (the  image  area  over  which  the  wavefront  error 
due  to  atmospheric  turbulence  is  correlated)  was  much  smaller  than 
the  actual  field  of  view  of  the  wavefront  sensor,  resulting  in  an 
averaging  of  the  measured  wavefront  error.  Under  these  conditions, 
the  "measured"  wavefront  error  did  not  represent  the  "true"  wavefront 
error  in  any  part  of  the  field  of  view  so  that  compensation  was  not 
effective. 

Analysis  of  anisoplanatic  effects  has  shown  that  with  normal  daytime 
turbulence,  the  isoplanatic  patch  near  the  zenith  may  only  be  1  to  2 
arc  seconds  in  radius,  rapidly  decreasing  as  the  eight  fifths  power 
of  the  cosine  of  the  zenith  angle. 
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This  small  isoplanatic  angle  produces  two  practical  difficulties: 


1)  If  the  wavefront  sensor  field  of  view  is  larger  than  the 
isoplanatic  patch,  the  degree  of  improvement  possible  in  image 
quality  is  very  limited. 

2)  Even  if  image  compensation  is  achieved,  it  will  be  effective 
over  a  very  small  area.  With  a  small  aperture  telescope  such  as 
the  30  cm  compensated  aperture  achieved  with  the  RTAC,  the  compen¬ 
sated  area  includes  only  a  small  number  of  resolution  elements. 

3.2  Second  Observing  Period 

3.2.1  Solar  Images 

The  RTAC  wavefront  sensor  as  configured  for  the  June  1979 
tests  was  capable  of  using  areas  of  the  solar  surface  of  15  x  15 
arc  seconds  down  to  5  x  5  arc  seconds  as  the  reference  source. 

To  improve  performance  with  small  isoplanatic  patches,  the 
RTAC  was  reconfigured  at  Itek  in  late  1979  to  reduce  the  wavefront 
sensor  field  of  view  to  about  2x2  arc  seconds.  A  second  series  of 
tests  was  conducted  using  the  76  cm  Vacuum  Tower  Telescope  at 
Sacramento  Peak  between  21  March  and  7  April  1980. 

This  modification  was  partly  successful  in  that  a  noticeable 
improvement  in  solar  image  quality  was  obtained  on  two  occasions, 
one  example  being  shown  in  Figure  5.  The  daytime  seeing  during  this 
period  varied  from  poor  to  moderate  and  it  was  necessary  to  use  a 
small  telescope  aperture  (0.3  m) .  Because  of  this,  the  smallest 
field  stop  size  usable  was  2.3  x  3.6  arc  seconds.  It  was  again 
found  that  on  most  occasions  the  RTAC  servo  loop  could  be  closed, 
reducing  the  indicated  electrical  error  signals  to  near  zero.  A 
corresponding  improvement  in  image  quality  was  seen  only  occasionally, 
and  could  not  be  predicted. 

3.2.2  Stellar  Images 

An  objective  of  the  second  series  of  tests  was  to  obtain 
isoplanatic  data  from  double  stars  by  using  one  star  as  the  reference 
and  measuring  the  image  quality  of  the  second  star  as  a  function  of 
the  angular  separation. 

To  facilitate  the  star  tests,  it  was  intended  to  connect  the 
RTAC  tip/tilt  outputs  to  the  tracking  system  of  the  76-cm  Tower 
Telescope,  which  normally  cannot  track  stars  accurately.  However, 
this  arrangement  was  not  successful  due  to  problems  in  one  of  the 
telescope  synchros  which  could  not  be  fixed  in  time.  We  therefore 
reverted  to  using  the  old  quad  PMT  star  tracker,  which  only  operates 
when  the  sky  is  completely  dark. 
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■4 


(a)  No  compensation 


Fig.  5  —  Real-time  compensation  of  sunspot  using  21-aUuator  deformable 
mirror  (30-cm  telescope  aperture,  18  x  13 -arc-sec  field  of  view,  1/15-sec 
exposure) 


We  also  found  that  the  RTAC  would  not  lock  on  with  3rd  Magnitude 
stars  due  to  the  optical  transmission  of  the  telescope  being  less 
than  anticipated.  These  two  factors  ruled  out  observation  of  double 
stars  from  night  into  daylight,  because  the  morning  stars  were  all 
3rd  magnitude.  Consequently,  our  observations  of  double  stars  were 
restricted  to  those  of  2nd  magnitude  visible  at  night,  which  were: 

a  UMa  ,  1  arc  second  separation.  Magnitudes  1.88,  4.82 
a  Gem  ,  2.2  arc  second  separation,  "  2.0,  2.8 

C  Ori  ,  3  "  "  "  "  1.91,  4.04 

y  Leo  ,  4.3  "  "  "  1.1  ,  3.4 

Images  of  several  single  bright  stars  such  as  Sirius,  Arcturus, 
and  Vega  were  recorded  in  addition  to  the  double  stars.  The  improve¬ 
ment  in  image  quality,  when  the  MPM  was  switched  on,  was  considerable, 
as  shown  in  the  intensity  profiles  of  Figure  6.  These  profiles  are 
based  on  the  integration  of  16  video  frames,  equivalent  to  about  1/4 
second  exposure  time. 

The  double  star  images  were  obtained  with  the  RTAC  operating  near 
its  photon  limit.  Consequently,  the  signal  to  noise  ratio  is  low 
and  the  compensation  not  as  effective  as  with  brighter  stars.  Inten¬ 
sity  profiles  of  a  Gem  (Castor)  and  y  Leo  (Algieba)  are  shown  in  Figures 
7  and  8.  In  the  case  of  Castor,  both  star  images  appear  to  be  equally 
changed  by  operation  of  the  MPM,  indicating  that  the  isoplanatic  patch 
exceeds  2  arc  seconds.  The  profiles  of  Algieba  clearly  show  sharpening 
of  the  brighter  star  (used  as  reference)  with  degradation  of  the  com¬ 
parisons  indicating  that  in  this  instance,  the  isoplanatic  patch  is 
less  than  4  arc  seconds  in  radius. 

3.3  Data  Summary 

An  important  output  of  the  Solar  Imaging  Experiment  is  the  large 
data  base  that  was  collected  by  recording  the  real-time  signal  outputs 
of  the  RTAC  when  observing  both  during  the  day  and  at  night.  These  data 
are  summarized  in  Table  3,  which  lists  both  the  videotape  cassette  and 
counter  numbers,  and  also  the  corresponding  digital  data  tape  sections. 
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Figure  6.  Intensity  Profiles  of  Arcturus 
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INTENSITY  PROFILES  OF  CASTOR 


56°  ELEVATION 
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Figure  7.  Intensity  Profiles  of  Castor 


INTENSITY  PROFILES  OF  ALGIEBA 
ALGIEBA 


Figure  8.  Intensity  Profiles  of  Algieba 


TABLE  3 

SAC  PEAK  OBSERVATIONS  MAR/APR  1980 
DATA  SUMMARY 


OBJECT 


Sunspot 

Sunspot 

Sunspot 

Sirius 

Castor 

Algieba 

Sirius 


Sunspot 


Arcturus 

ft 

Sirius 


Zeta  OR! 

Castor 

Algieba 


DATE/TIME 

HRS. 


3/28  10:00 

3/28  11:30 

3/29  16:15 

3/30  20:00 

3/30  22:30 

3/30  23:00 


VIDEOTAPE 


DATA  TAPE 
# 


3/31  21:30 

3/31  21:30 

3/31  21:30 

3/31  21:30 

3/31  21:30 

3/31  21:30 

4/03  11:30 

4/94  07:30 

4/04  09:30 

4/04  11:00 

4/05  09:30 

4/06  02:45 

4/06  03:00 

4/06  19:45 

4/06  20:00 

4/06  21:00 

4/06  21:30 

4/06  22:30 


156-84 

2 

120-140 

3 

145-166 

4 

183- 188 

5 

188-218 

5 

234-280 

5 

000-015 

6 

- 

6 

- 

6 

- 

6 

- 

6 

- 

6 

032-097 

- 

097-117 

- 

125-134 

7 

1  4 1  - 1  5  3 

/ 

200-210 

- 

032-050 

8 

050-080 

- 

083-125 

9  0  to 

1  minll 

135- 148 

- 

207-220 

9  1  min 

10  to  3  min' 

244-260 

9  3  min 

40  to  6min  40 

320-330 

- 

353-368 

REMARKS 

mi 

Best 

^Sunspot  images 
turbulence 

Best 

turbulence 

Tilt 

only 

Tilt 

only 

TABLK  3 


(Continued) 


OBJECT 


DATE/TIME 

HRS. 


5  090-119 


120-133 
202-  250 


250-265 

317-328 

328-352 


205-  214 
231-250 
250-266 


END  OF  OBSERVATION 


VIDEOTAPE  DATA  TAPE 

#  INDEX  # 


5  015-031 


REMARKS 

015-020  Good  031 

No  turbulence,  RTAC 
Baseline 

Best  sections  90-91,9 
113-114,  117-119  (BWl 

Best  224-226 
56°  EL  236-250 

BWl  -  Best  ACC  OFF 
BW2  -  Not  as  good 

210-212  OK  30°  El<  v 

20°  Elev 


Software  to  process  these  data  tapes  has  been  generated  on  the 
contract.  The  specific  programs  developed  are: 

a)  Summary  plot  of  4  channels  of  data  from  the  entire  data 
tape,  with  Universal  time  scale. 

b)  Power  Spectrum  computation  and  plot  for  designated  2- 
second  periods  of  any  channel. 

The  following  specific  items  have  been  delivered  under  the 
contract : 

1.  Documentation  of  the  Power  Spectrum  Software. 

2.  Sample  test  run  producing  both  universal  time  plots  and 
power  spectrum  plots. 

3.  Source  code  for  Power  Spectrum  Software  excluding  the 
FILENO,  FOURT,  and  READIN  routines.  Both  FILENO  and  FOURT  may  be 
obtained  from  the  SHARE  Fortran  Project  and  READIN  was  developed 
by  ITEK  and  cannot  be  distributed. 

4.  Documentation  for  FILENO,  FOURT,  and  READIN. 

5.  Magnetic  Tape  containing  source  code  for  Power  Spectrum. 

6.  Magnetic  Tape  containing  power  spectrum  data. 

Examples  of  the  data  recorded  are  shown  in  Figure  9  and  a  typical 
Power  Spectrum  plot  is  shown  in  Figure  10. 

Results  of  the  data  analysis  performed  to  date  are  shown  in 
Table  4,  Figure  11  and  Figure  12.  Table  4  lists  the  temporal  band¬ 
width  of  the  drive  signals  for  tip-tilt  and  MPM  for  both  daylight  and 
nighttime  conditions,  together  with  the  root  mean  square  values  of 
the  drive  signals.  Also  listed  is  a  single  measurement  of  the  band¬ 
width  of  the  intensity  (scintillation)  signal  for  a  star  (Sirius)  at 
night . 

The  tip  tilt  data  show  little  difference  between  day  and  night 
operation,  the  bandwidth  measured  at  101,  power  ranging  from  about 
10  to  40  Hz. 

A  much  greater  difference  is  seen  in  the  MPM  drive  signal  power 
spectra  between  day  and  night  oneration,  the  daytime  10%  power  band¬ 
width  ranging  from  60  to  120  Hz,  while  at  night  it  lies  between 
20  and  55  Hz.  The  bandwidth  is  of  course,  highly  dependent  on  wind 
velocity. 
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FIGURE  9  REAL-TIME  DATA  PRINT-OUT,  CLOSED-LOOP  OPERATION. 

ITEK  DATA  TAPE  NO.  3A,  SECOND  2-SECONDS  OF  DATA. 
SUNSPOT,  28  MARCH,  1980,  1130  HRS.  LOCAL  TIME. 
COMPENSATED  APERTURE  30  CM  DIAMETER. 


POWER  NCHRN  =  11 


igure  10  -  Power  Spectrum  of  Actuator  Drive  Signal 


23 


TABLF.  4 


DATA  ANALYSIS  SUMMARY 
Temporal  Bandwidth  and  Drive  Distances 


A.  TILT 


Date 

and 

Time  Conditions 

1 

!Tap. 
I  ff 

5/28/80 

1130 

SUNSPOT  58° 

5 

i 

,  3A 

4/06/80 
02  30 

ARCTURUS  70° 

,  8 

4/06/80 

2000 

SIRIUS 

30° 

i 

,  9A 

B.  MPM 

3/ 28/ 80| 
1130 

i 

SUNSPOT 

58° 

!  3A 
;  2 

3/31/80! 

2120 

SIRIUS 

20° 

j  6A 

4/06/80 

2000 

SIRIUS 

30° 

9  A 

1 

3/31/80 

2130 

SIRIUS 

18° 

| 

;  6B 

3/31/80 

2  140 

SIRIUS 

17° 

3/31/80 
2150  i 

SIRIUS 

16° 

l6D 

(’.  I NTFNS  1  TY 

4/06/80 

2000 

SIRIUS 

1  9A 

1 

1  ... 

I  Aper- 

Bandwidth _ j  ture 

Cutof f _ 10i~j  Meters 


X 

50 

31 

.  0.3 

Y 

60 

31 

X 

100 

40 

0.4 

Y 

55 

20 

1 

X 

35 

12 

•  0.46 

Y 

18 

8 

168 

124 

0.3 

123 

62 

120 

50 

0.3 

125 

55 

,0.46 

100 

50 

0.35 

100 

45 

0.41 

90 

20 

j  0.53 

1 

104 

,  34 

1 

1 

+  1 
+  0 

+  2 
+  1 


+  0 


+  n 


+  o 


+  o 


Drive  Angle 
or  Distance 


78  arc  sec  RMS 
87  arc  sec  RMS 

16  arc  sec  RMS 
71  arc  sec  RMS 


48uM  RMS 

48pM  RMS 
70UM  RMS 

50uM  RMS 
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Correlation 


Meters 


Spatial  Correlation,  Night  Observations 
from  Tape  No.  6,  Sirius,  31  March  1980. 

Each  point  is  a  2-second  average. 

A  *  20°  Elevation  B  =  18°  Elevation 
C  *  17  Elevation  D  =  16°  Elevation 


Fig.  11 


Spacing,  Meters 


Intensity  Correlation,  Night  Observations 
from  Tape  No.  9A,  Sirius,  6  April  1980. 

Each  point  is  a  2-second  average. 
Elevation  angle  30  . 


The  spatial  correlation  of  the  MPM  drive  signals  is  shown  in 
Figure  11,  for  night-time  operation  using  2-second  samples,  taken 
over  a  30  minute  period.  The  distance  for  50%  correlation  varies 
between  5  and  15  cm.  Random  changes  between  samples  appear  to 
dominate  any  systematic  change  with  elevation  or  time.  Spatial 
correlation  of  the  photometric  intensity  of  the  received  wavefront 
is  shown  in  Figure  12.  The  50%  correlation  distance  was  between 
11  and  15  cm,  somewhat  greater  than  that  for  the  wavefront  phase. 
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4. 


ANALYSIS 


The  optical  effects  of  atmospheric  turbulence  can  be  characterized 
by  three  parameters:  the  spatial  coherence  length,  rQ;  the  iso- 
planatic  angle  0Q ;  and  the  temporal  bandwidth  fo.  The  temporal 
bandwidth  is  mainly  determined  by  wind  velocities  and  is  normally 
in  the  order  of  a  few  hundred  Hertz.  To  keep  the  compensation 
error  due  to  temporal  bandwidth  small,  the  wavefront  measurement 
time  and  the  servo  response  time  must  not  exceed  a  few  millisecond' 
With  a  bright  object  such  as  the  sun,  this  presents  no  problem 
because  the  large  photon  flux  allows  the  wavefront  measurement  to 
be  made  with  a  high  signal  to  noise  ratio  using  an  integration  time 
of  less  than  a  millisecond.  Also,  the  response  of  current  deform¬ 
able  mirrors  is  well  above  1  KHz.  The  atmospheric  bandwidth,  fo, 
is  therefore  not  a  practical  limitation,  so  we  will  concentrate 
attention  on  the  remaining  two  parameters. 

4.1  Coherence  Length 

The  usually  accepted  definition  of  the  turbulence  coherence 
length,  r  ,  due  to  Fried,  is  in  the  form  of  the  turbulence  structure 
funct ion ,  ° 

<  ( 0 1  -  0  2 )  >  =  D(r)  =  6.88  (— !-)  /  radians 

ro 

which  defines  the  mean  square  phase  variance  between  two  points  on 
the  wavefront  separated  by  distance  r.  The  mean  square  phase 
variance  is  a  measure  of  the  total  integrated  turbulence  of  the 
atmosphere  along  one  line  of  sight.  Due  to  the  exponential  fall- 
off  of  turbulence  strength  with  altitude,  the  lower  layers  of  the 
atmosphere  have  a  dominent  effect  on  r  .  Because  of  solar  heating 
near  the  ground,  a  strong  diurnal  effe8t  is  present.  This  effect 
is  shown  in  Figure  13  which  is  derived  from  data  obtained  by 
Walters  et  al .  at  White  Sands  (Reference  2).  Note  that  relatively 
high  values  of  r  (low  turbulence)  occur  during  the  two  periods  of 
temperature  equilibrim  around  sunrise  and  sunset. 

In  addition  to  the  diurnal  variation,  the  effective  value  of  r  , 
normally  specified  at  zenith,  is  a  function  of  zenith  angle,  t? 

,  .  3/5 

r  (  Q  ~  r  cos  '  c, 
o  o 

at  low  elevations,  therefore,  the  effective  value  of  r  may  be 
very  small  (high  turbulence)  even  though  the  zenith  seeing  may 
be  good. 

'! o  determine  the  effective  solar  seeing,  we  must  therefore  include 
both  diurnal  and  zenith  angle  effects,  as  shown  in  Figure  14.  The 
best  periods  of  seeing  occur  in  early  morning  and  late  afternoon, 
but  it  should  be  noted  that  the  values  of  r  obtained  (5  to  6  cm  ) 
even  then  are  only  hali  of  those  encountereS  at  night  at  good 
astronomical  sites. 
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0  2  4  6  8  10  12  14  16  18  20  22  24 

Time  of  day 

Fig.  13  Diurnal  variation  of  turbulence  coherence 
length  at  sites  A,  3,  and  C 
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Average  of  three  sites 
Best  of  three  sites 


4  6  8  10  12  14  16  18  20 

Time  of  day 

Fig.  14  Variation  of  effective  solar  coherence  length 
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to  a  wave- 
wavefront 


where  d  =  spacing  of  deformable  mirror  actuators 

K  -  fitting  coefficient  of  mirror 

The  factor  K  depends  on  the  shape  or  influence  function  of  the 
actuators;  for  Itek's  monolithic  piezoelectric  mirror  it  has  the 
value  of  0.01  waves.  The  value  of  o ~p  can  always  be  reduced  by 
decreasing  d.  In  practical  systems,  the  number  of  actuators  is 
fixed  and  d  can  only  be  reduced  by  using  a  smaller  aperture. 
However,  this  limits  the  angular  resolving  power  of  the  telescope 
so  a  compromise  is  necessary.  For  the  experiments  described  in 
Section  3,  an  optical  zoom  system  was  used  between  the  telescope 
and  the  RTAC  to  vary  the  effective  aperture  between  30  cm  and 
60  cm,  corresponding  to  values  of  d  between  6  cm  and  12  cm,  in 
order  to  match  the  expected  range  of  rQ. 

4.2  Isoplanatic  Patch 

rf  we  consider  a  layer  of  turbulence  with  a  specific 
coherence  length,  then  it  is  obvious  that  the  isoplanatic  angle 
over  which  the  wavefront  is  correlated  will  be  largest  when  the 
layer  is  near  the  telescope  aperture  and  smallest  when  it  is  far 
from  the  aperture.  The  isoplanatic  angle,  0  is  therefore 
determined  by  the  vertical  distribution,  or  profile,  of  atmos¬ 
pheric  turbulence,  with  the  higher  layers  having  the  dominant 
effect.  The  isoplanatic  angle  is  largely  independent  of  r  ,  and 
may  be  relatively  independent  of  diurnal  effects.  The  isoplanatic 
angle  f  Q  can  he  defined  in  a  similar  way  to  coherence  length: 

0  5/3  2 

D(f  j  =  6.88  f — -  )  radians 

O 

which  defines  the  mean  square  phase  variance  between  two  rays  at 
an  angle  <>.  Measurements  made  by  various  workers  (References  3  and  4) 
suggest  that  the  value  of  o  as  defined  above  is  generally  in  the  region  of 
3  to  6  arc  seconds,  measurea  at  the  zenith. 

The  isoplanatic  angle  is  strongly  dependent  on  zenith  angle,  the  relation 
being: 

8 /S 

no  ^  =  eo  cos  r‘ 


The  practical  importance  of  the  coherence  length  r 
front  compensation  system  is  that  it  determines  th§ 
fitting  error,  which  is  given  by: 

2  a  5/3 

o  P  =  K  ( — )  ' 

F  r 


The  isoplanatic  angle  o  defined  above  is  a  property  of  the  atmospheric 
turbulence  profile,  independent  of  any  compensation  system  and  pertains  to 
single  rays  traversing  the  atmosphere. 


When  a  point  source  is  used  for  wavefront  sensing,  the  measured 
wavefront  is  exactly  correct  in  the  direction  of  the  reference 
source.  A  wavefront  error  variance  proportional  to  o5'  is 
incurred  as  the  field  angle  0  increases. 

With  an  extended  object  such  as  the  sun,  a  complication  arises 
in  that  a  finite  area  of  the  object,  encompassing  some  surface 
detail  such  as  granulation,  must  be  used  as  the  reference  source. 
The  wavefront  is  then  compensated  with  regard  to  the  average 
error  over  this  finite  reference  area.  It  is  therefore  of 
interest  to  determine  the  effect  of  the  size  of  this  reference 
area  on  the  anisoplanat ic  error  of  the  compensated  image. 


The  results  of  an  analysis  made  at  Itek  are  summarized  in  Figure 
IS  which  shows  the  mean  square  wavefront  error  a2  5due  to  aniso- 
planatism,  normalized  by  the  function  6.88  (0p/oo)  /  ,  as  a 
function  of  D/Dr,  where  0  is  the  field  angle  and  0p  is  the  radius 
of  the  area  usea  for  wavefront  reference.  The  main  effect  of  using 
a  finite  area  rather  than  a  point  source  for  wavefront  reference 
is  to  increase  the  wavefront  error  in  the  center  of  the  field-of- 
view  and  to  reduce  it  at  the  edge  of  the  field. 


It  should  be  noted  that  this  effect  is  a  property  of  the  atmos¬ 
phere  and  is  present  even  with  a  "perfect”  wavefront  sensor. 

Also,  the  minimum  size  of  the  wavefront  reference  field  is  deter¬ 
mined  by  the  scale  of  the  solar  detail  (such  as  the  granulation 
pattern);  it  is  not  a  discretionary  variable  and  is  the  same  for 
any  type  of  wavefront  sensor. 


To  determine  the  overall  image  resolution  as  a  function  of  field 
angle,  both  anisoplanat  ic  and  fitting  error  effects  must  be  taken 
into  account.  Figure  16  shows  plots  of  on-axis  image  intensity 
versus  field  angle  for  typical  conditions  when  viewing  the  sun 
near  the  zenith  and  at  30  degrees  elevation.  At  the  zenith,  good 
compensation  is  obtainable  over  a  field  of  about  2  arcseconds  in 
diameter.  At  30  degrees  elevation  however,  the  isoplanatic  angle 
has  shrunk  to  one  third  of  its  zenith  value  and  r  to  two  thirds 
of  its  zenith  value,  resulting  in  a  drastically  reduced  compensa¬ 
tion  patch  less  than  1  arcsecond  in  diameter. 


It  is  of  interest  to  determine  how  the  size  of  the  compensated 
solar  image  varies  during  the  day.  The  compensated  image  size  is 
defined  arbitrarily  as  the  radius  at  which  the  wavefront  error  is 
0.20  waves  rms ,  giving  an  on-axis  image  intensity  of  0.2  of  its 
unaberrated  value.  Using  the  diurnal  variation  of  rQ  from  Figure 
14  the  results  shown  in  Figure  17  are  obtained. 


These  results  show  that  with  solar  image  compensation,  the  largest 
compensated  area  is  obtained  around  midday.  During  the  periods  of 
good  uncompensated  seeing  in  early  morning  and  late  afternoon,  the 
size  of  the  compensated  patch  is  too  small  to  be  useful. 
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This  analysis  provides  an  insight  into  the  results  obtained 
during  the  solar  image  compensation  experiment,  as  discussed 
in  the  following  section. 


5.  CONCLUSIONS 


When  observing  stars,  the  real-time  compensation  system  operated 
as  expected,  sharpening  the  image  and  producing  a  significant  in¬ 
crease  in  the  peak  intensity.  The  quality  of  the  compensated 
image  was  generally  consistent  with  the  fitting  error  limitations 
due  to  the  relation  of  actuator  spacing,  d  to  atmospheric  r  ,  as 
described  in  Section  4.1.  Observations  of  double  stars  were  some¬ 
what  restricted  due  to  magnitude  limitations.  The  data  obtained 
indicate  that  the  isoplanatic  patch  radius  at  night  was  between 
2  and  4  arc  seconds  measured  near  the  zenith. 

The  results  obtained  when  observing  the  sun  were  not  as  expected. 

On  most  occasions  the  RTAC  appeared  to  operate  normally;  the  active 
mirror  feedback  loops  were  converged  and  stable;  the  indicated  re¬ 
sidual  error  was  reduced  to  a  small  value;  the  real-time  data  outputs 
from  the  RTAC  (such  as  the  MPM  drive  signals  shown  in  Figure  9)  were 
exactly  as  expected  in  terms  of  drive  distance  and  power  spectrum. 
However,  a  significant  increase  in  image  resolution  was  obtained  only 
on  rare  and  unpredictable  occasions. 

While  initially  these  results  were  puzzling,  we  believe  that  they 
can  be  explained  in  terms  of  the  system  parameters,  particularly 
the  values  of  r  and  0  in  relation  to  the  aperture  size  and  number 
of  actuators  onthe  RTAC. 

figure  17  shows  that  under  average  conditions  (r  =  4  cm,  0*3 
arc  seconds'!  with  an  actuator  spacing  of  6  cm  ana  a  wavefront  sensor 
field  stop  radius  of  1.5  arc  seconds,  the  compensated  image  radius 
(for  0.2  waves  rms  error)  is  less  than  1  arc  second.  During  the 
periods  of  good  solar  seeing  around  8  to  9  AM,  the  compensated  image 
radius  is  less  than  0.5  arc  seconds. 

With  a  21-elemcnt  active  mirror  having  5  actuators  across  the  dia¬ 
meter,  the  6  cm  actuator  spacing  corresponds  to  a  compensated  tele¬ 
scope  aperture  of  30  cm;  the  classic  resolving  power  of  this  aperture 
at  X  =  500  nm  is  0.4  arc  seconds.  The  number  of  resolution  elements 
across  the  diameter  of  the  compensated  field  is  therefore  between  2.5 
and  5.  Such  a  small  compensated  area  is  difficult  to  distinguish  vis¬ 
ually  in  a  video  image  frame  covering  15  arc  seconds,  especially  when 
the  image  is  of  low  contrast. 

On  the  occasions  when  a  noticeable  improvement  in  image  quality  was 
obtained  (figure  5,  for  example),  the  area  of  improvement  was  relatively 
large  (about  6  arc  seconds  in  this  case)  indicating  that  on  this  parti¬ 
cular  occasion  0Q  was  much  larger  than  its  average  value. 

The  results  of  the  experiment  are  generally  confirmed  by  the  analysis, 
which  shows  that  even  with  a  perfectly  functioning  real-time  compensation 
system,  the  size  of  the  compensated  area  on  the  sun  is  determined  pri¬ 
marily  by  atmospheric  conditions  and  is  normally  only  one  or  two  arc 
seconds  in  diameter. 

The  benefits  of  real-time  compensation  become  more  apparent  as  the 
telescope  aperture  is  increased,  because  a  larger  number  of  image  ele¬ 
ments  can  then  be  resolved  within  this  compensated  area. 
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